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A new semiclassical approach to linear (L) and nonlinear (NL) one-dimensional Schrodinger 
equation (SE) is presented. Unlike the usual WKB solution, our solution does not diverge at the 
classical turning point. For LSE, our zeroth-order solution, when expanded in powers of h, agrees 
with the usual WKB solution e^^^^ up to 0{h'^), where 5* and v are the classical action and 
the velocity, respectively. For NLSE, our zeroth-order solution includes quantum corrections to the 
, Thomas-Fermi solution, thereby giving a smoothly decaying wave function into the forbidden region. 

^ ■ 03.65.Ge, 03.65.Sq 

■ The WKB method allows us to derive expressions for various quantum-mechanical quantities when the action is 
^ much larger than fi, and has been widely used in many subfields of physics and chemistry [Q-^ . Nevertheless, the WKB 

method has a serious flaw of divergence at the classical turning point, which originates from the fact that the classical 
trajectry is taken as the zeroth-order solution. Considerable efforts have been devoted to overcome this problem, 
^sj among other things, the complex method [^j^ and the uniform approximation |^,^. However, the former method 
04 cannot describe the wave function at the turning point, while the latter one requires us to resort to a comparison 
function |0] that mimics the original potential and to which an analytic solution is available. This Letter presents a 
^ ■ new semiclassical method whose zeroth-order solution is constructed upon a trajectry with quantum corrections such 
. that the aforementioned divergence is absent. Moreover, our method gives a solution to the nonlinear Schrodinger 
' equation (NLSE) on an equal footing. 

' We begin by reviewing the WKB method for the one-dimensional linear Schrodinger equation (LSE) — ^ = 
I , Etp. Rescaling the length and the energy in units of / and /2'mP, respectively, where I is a characteristic length 

■ scale of the potential V{x), LSE takes the form 

o : 
' 



O 



■4)" + Vi; = Eil; , (1) 



^ , where the prime denotes the differentiation with respect to x. In Eq. (|1|), the length x and the energies E,V are 
I ' proportional to hP and h~'^ , respectively. We introduce the action function ip{x) through ^{x) = e'^'^^^ , where the 
^ , action is measured in units of %. In terms of ^{x), Eq. (Q) takes the form 



i^r + iE-V)^-^". (2) 
The zeroth-order WKB solution (^wkb o > which is obtained by neglecting ip" in Eq. (||), satisfies the following equation: 

(^Wkb.o)' + {E-V) = 0. (3) 

The classical turning point xt'^) is determined by ^(x*^'^^) = E. Incorporating the effect of ip" perturbatively, the WKB 
solution <(5^j^3 takes familiar forms as 

We note that the WKB solution diverges at a;*^*^^. This divergence originates from the fact that the classical trajectry 
is taken as the zeroth-order solution. Our strategy is to incorporate quantum corrections in the zeroth-order solution 
so as to remove the divergence. 

To carry out this program, we differentiate both sides of Eq. (^) with respect to x, having 

2^V'-r = V, (5) 
and substitute tp" in Eq. (H) into Eq. (|), obtaining 

(<pT + (£^- W + = ■ (6) 
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The Schrodinger equation (|^) is sufRcient for Eq. to hold, but it is not necessary. To show this, we add </?V" to 
both sides of Eq. and integrate the resulting equation with respect to x. We then get 



ge'^ - {ip'f + {E-V) + if" 
where g \s a, constant of integration. Equation can be rewriten as 

- ip" + FV' + g^^ = EiP . 



(7) 



(8) 



This is nothing but NLSE which includes LSE as a particular case of g = 0. Throughout this Letter, we assume that 
for NLSE the wave function is normalized to unity dx \ij}{x)\'^ = 1, while for LSE the normalization constant is 
left arbitrary. 

Our zeroth-order solution (^q, which is obtained by neglecting ip'" in Eq. (p), satisfies the following cubic equation 



{f'^f + Zpf>'^ + 2g = , 



(9) 



where p = [E — y)/3 and q = V /4. Comparing Eqs. (||) and (||) with Eqs. (^) and (^), respectively, we see our 
method is a natural extension of the WKB method. Equation is the fundamental equation of our new semiclassical 
method. As we shall see later, the solutions of Eq. (^) lead to divergence-free wave functions. 

A root of the determinant of Eq. D = + =0, determines our turning point x — x^'^\ which includes 
quantum corrections. In fact, the difference between x'--'^'^ and a;*-'^-' may be expanded as 



,(q) 



(c) _ 
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where 



~ [2y'(a;W)]4/3 



and S 



We first discuss the zeroth-order solution. In the allowed region, the cubic equation has the following three 
solutions: 



^K + ik = if'^ 



where 




(11) 



k{x) = 

k{x)= ^ + V:D)1/3_ (q_7;D)l/3' 



(12) 



When expanded in powers of h, Eq. (11) reduces to 
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where —k and 93 respectively, correspond to the Thomas-Fermi solution 
0{h°). Substituting Eq. ^ into Eq. (0), we obtain 

( E^V+ 0{h°) for ip'^ = -K , 

[ 0{h°) for ifi'^ ^ ip'_ . 

From Eq. (|l^), it is clear that —k and Lp'j^ give a solution to NLSE and those to LSE, respectively. 
In the forbidden region, the cubic equation ^ has the following three real solutions : 



(13) 



up to 



(14) 
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-K = ^2y/—pcos ( iArctan- 



^'a={ X+ = T2\/^cos(iArctan^ + f ) , (15) 



= =f2\/^cos ( iArcta 



-D 2tt ' 
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where the — and + signs correspond to V > and V < 0, respectively. When expanded in powers of h, Eq. ([l5| ) 
reduces to 

(p;,^ <{ + j^ + oih')= ±Vv-E - i^ + oih') , (16) 

2q _ 
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where —k and x+ correspond to the WKB solutions (Q) up to 0(?i"). Substituting Eq. (|T^) into Eq. (|^), we obtain 

r 0(^°) for ip'o = -K , 

9^^ = \ 0{h") for ^'o^x'+, (17) 

[ + 0{h°) for (^[, -x'_ . 

To proceed further with our analysis, we shall assume that V'{x) > (Vx). Then the wave function must decay 
to zero as a; ^ oo, i.e., lim ■(/'(x) = 0. Consequently, in the forbidden region, —k must be chosen as the zeroth-order 

X — *OD 

solution (p'q for both cases of LSE and NLSE. Note that, for NLSE, gip'^ = 0{lf) in Eq. (0) does not mean that g = 
but that the wave function is well damped in the forbidden region. 

We are now in a position to construct our zeroth-order wave functions to LSE and NLSE. For NLSE, the zeroth-order 
wave function is 

^o{x) ^ N exp ( ~ [ dx' k{x')] , (18) 



where TV is a normalization constant. Note that this single solution covers both allowed and forbidden regions. 
For LSE, the zeroth-order wave function can be described, in general, as follows 

_ r exp (/;,,, dx' ip'+{x')) + A_ exp (/;,,, dx' ip'_{x')) ^ M^) < ^('^)) , 

^''^'^^ \ exp (- dx' k{x')) EE tPuix) (x > x(^^) . 

To determine the relation between constants and A_, we note that, near x^'^\ LSE (^ reduces to 



(19) 



dSp 



TJTT + = ' (20) 



where ^ = [y'(a;('=))]i/3( X — x^'^'). The exact solution of Eq. (GO) that satisfies the boundary condition lim tp{^) — 

— ^^oo 

is the Airy function Ai(^) [|l^. In the region |^| 3> 1, the asymptotic forms of Ai(^) are 

f;fe(-er^/^-(|(-o;/^-,f) (e<o), 

l2^r^/^ exp(-|C3/2) (^>0). ^^^^ 

On the other hand, for the linear potential, where p — — JV'(x('^))]^/'^^/3 and q = y(a;('^^)/4, one can integrate ip'j. 
and K in Eq. ([l^); then the asymptotic behavior of Eq. ( [l^ ) for |^| ^ 1 becomes 

r Vi - (2e)-2/32V3(_e)-i/4 [A^ exp {~^ (|(-0^/^ + f )) + ^- exp (§(-0^/^ + f ))] 
^0 = <^ (C < 0) , (22) 

i 1^11 ^ i(2e)5/62i/3^-i/4 exp (-f^^/^) (C > 0) . 

For Eq. ( p2| ) to match Eq. (|2l]), we must choose = — 1. The zeroth-order solution ( p^ can thus be described 

as 
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i/,o(a;) ^ I V'l = ^exp (i /^"i,, dx' k{x')) sin (^jf^ dx' fc(x')) {x < Xq) , 
[ = exp (- /^^(,) dx' k{x')) {x > xq) , 



(23) 



where A = ^2iA±. For and -i/jn and their first derivatives to be continuous at a certain point xq (< x^'^^), i.e. 



2k(xo) 
3k.(xo) 



and 



A 



"01 (2^0) = "011(2^0) a-iid ^{(xo) = ^^(xo), xo and A must satisfy tan ( J^^ dx k{x)j 

I J^^ dx k{x) j cosec ( J^^ dx k(x)j . These simulutaneous equations can be solved approximately, giving xq 

xi'^') + 0{e) and A ~ (2e)^^^ + 0(e). For the linear potential (e = 0), the asymptotic behavior of ■01 (^ii) in Eq. (|22 
for ^ — > —00 (+00) with A± = ^{2e)^^^ /2i exactly agrees with that of Ai(^) in Eq. (21) except for an overall factor. 

Our zeroth-order wave function ( p3| ) vanishes at x^'^' for any potential, which gives us a new quantization condition 
(compared with the WKB quantization condition) 



^§dx k{x) = i /""(^j dx k{x) =n + l 

(Our method), 
{x) = - /_^(^) dx fcwKB(a;) =n + ^ 

(the WKB method). 



/ ftwKB 



(24) 



where k{x) is given in Eq. (|T^), fcwKB(2;) = \/E — V{x), n = 0, 1, 2, • • • is the number of nodes of the wave function, 
and the suffixes L and R denote the left- and right-side turning points, respectively. We note that k(x) and x^'^-' 
include quantum corretions, which accounts for an extra 1/2 in our quantization condition. 

We have so far ignored the term ^p'" /2 on the right-hand side of Eq. (||). The effect of this term may be evaluated 
perturbatively, giving [O 
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The first-order wave function '\\}\ (x) can be constructed by using kp'^ + ip\ instead of ip'^ . 

We finally apply the zeroth-order solution iiQ^x) and the first-order solution i\^\{x) to a linear potential (Figs. 1 and 
2), a parabolic potential V = with E = 17 (Figs. 3 and 4), and the Morse potential V ~ 900[exp(— 2a;) — 2 exp(— a;)] 
with E — —(39/2)^ (Fig. 5). We also compare them with the exact solution and the WKB solution or a solution in 
which the Thomas-Fermi and WKB ones are combined. 

As shown in the figures, the zeroth-order wave function ipQ does not diverge at the turning point but the error 
is discernible around the classical turning point x'^'^\ This small discrepancy is drastically improved by using the 
first-order wave function ipi. The error associated with the zeroth-order solution is estimated to be [ pT| 



0.028 - 0.36e + O(e2,(5) . 



(26) 



Equation (|2^) implies that our zeroth-order wave functions ( |1§| ) and (^3|) are suitable for e > rather than e < 0; 
when e < 0, the wave function is unstable against a small perturbation. 

In conclusion, we proposed a divergence-free semiclassical method that enables us to solve LSE and NLSE on an 
equall footing. Our zeroth-order solution is constructed upon a trajectry that includes quantum corrections and 
therefore allows a rapidly converging perturbative expansion of the wave function. Our method may therefore be 
applied to drastically improving related semiclassical methods such as instantons ||l^ and periodic orbital theory [p^ . 
The results of these subjects will be reported elsewhere. 

This work was supported by a Grant-in- Aid for Scientific Reseach (Grant No. 11216204) by the Ministry of Edu- 
cation, Science, Sports, and Culture of Japan, and by the Toray Science Foundation. 
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FIG. 1. Zeroth-order solution ipn (dashed curve) and first-order solution ipi (dot-dashed curve) to LSE for a linear potential, 
—d^ip/d^^ + ^Tp = 0. The exact solution (solid curve) and the usual WKB solution (dotted curve) are superimposed for 
comparison. The region around the turning point is enlarged in the inset. 
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FIG. 2. Zeroth-order solution ipo (dashed curve) and first-order solution tpi (dot-dashed curve) to NLSE for a linear potential, 
—(fip/d£,'^ + S,ip + = 0. The exact solution (solid curve) and a solution in which the Thomas- Fermi and WKB ones are 
combined (dotted curve) are superimposed for comparison. The region around the turning point is enlarged in the inset. 
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FIG. 3. Solutions to LSE for a parabolic potential, V = with E = 17 and e = 0.0476. The notations are the same as 
FIG. 1 . The region around the turning point is enlarged in the inset. 
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FIG. 4. Solutions to NLSE for a parabolic potential, V = x with E = 17 and e = 0.0476. The notations are the same as 
FIG. 2. The region around the turning point is enlarged in the inset. 
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FIG. 5. Solutions to LSE for the Morse potential, V = 900[exp(-2a;) - 2exp(-a;)] with E = -(39/2)^, el = 0.0982 and 
eR = —0.0394. The notations are the same as FIG. 1. The regions around the two turning points are enlarged in the insets. 
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